arXiv: 1505.06077vl [hep-ph] 22 May 2015 


BONN-TH-2015-07, CERN-PH-TH-2015-123 


Precise determination of the Higgs mass 
in supersymmetric models with vectorlike tops 
and the impact on naturalness in minimal GMSB 

Kilian NickeQ 

Bethe Center for Theoretical Physics & Physikalisches Institut der Universitdt Bonn, 

53115 Bonn, Germany 

Florian StautH] 

Theory Division, CERN, 1211 Geneva 23, Switzerland 

Abstract 

We present a precise analysis of the Higgs mass corrections stemming from vectorlike top partners 
in supersymmetric models. We reduce the theoretical uncertainty compared to previous studies in 
the following aspects: (i) including the one-loop threshold corrections to SM gauge and Yukawa 
couplings due to the presence of the new states to obtain the DR parameters entering all loop 
calculations, (ii) including the full momentum dependence at one-loop, and (iii) including all two- 
loop corrections but the ones involving and g 2 . We find that the additional threshold corrections 
are very important and can give the largest effect on the Higgs mass. However, we identify also 
parameter regions where the new two-loop effects can be more important than the ones of the 
MSSM and change the Higgs mass prediction by up to 10 GeV. This is for instance the case in 
the low tan/3, small Ma regime. We use these results to calculate the electroweak fine-tuning of 
an UV complete variant of this model. For this purpose, we add a complete 10 and 10 of SU{5) 
to the MSSM particle content. We embed this model in minimal Gauge Mediated Supersymmetry 
Breaking and calculate the electroweak fine-tuning with respect to all important parameters. It 
turns out that the limit on the gluino mass becomes more important for the fine-tuning than the 
Higgs mass measurements which is easily to satisfy in this setup. 
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I. INTRODUCTION 


The discovery of the Higgs boson with a mass of about 125 GeV mg has a strong 
impact on the parameter range of supersymmetric (SUSY) models, especially as its mass 
value is turning into a precision observable with an uncertainty below 1%. In particular, in 
constrained versions of the Minimal Supersymmetric Standard Model (MSSM) large regions 
of the parameter space are not consistent with this mass range |3]. This is in particular 
the case for models where SUSY breaking is assumed to be transmitted from the hidden to 
the visible sector via gauge interactions like in minimal Gauge Mediated SUSY Breaking 
(GMSB). Even relaxing the predictive boundary conditions of a constrained model and 
considering the phenomenological MSSM with many more parameters at the SUSY scale, it 
is still rather difficult to find regions with the correct Higgs mass. Either, a very large mixing 
in the stop sector or heavy stop masses are needed to push the Higgs mass to the desired 
range HHIH]. However, the large stop mixing with light stops turns out to be dangerous 
because of charge and colour breaking minima [T91I23] . On the other side, very heavy stops 
introduce again a hierarchy problem which SUSY was supposed to solve. The question about 
naturalness and fine-tuning is even more pronounced in regions the small tan (3 region which 
recently gained some interest because of Higgs fits dMZ!: in these regions the tree-level 
Higgs mass is suppressed by a factor cos(2/3) and even much bigger loop corrections are 
needed than for larger values of tan (3. 

A widely studied ansatz to solve this tension and to reduce the necessary fine-tuning in 
SUSY models is to enhance the Higgs mass already at tree-level. For this purpose models are 
considered which give new F- |2ni32] or H-term contributions to the Higgs mass |33ti39] . 
The fine-tuning in these models is often better by a few orders compared to the MSSM. 
Alternatively, one can also consider models which give new loop-corrections due to the 
presence of additional large couplings to push the Higgs mass. This happens for instance 
in inverse-seesaw models HOI HI] or models with vector-like quarks j32HS3] at the one-loop 
level, or in models with trilinear i?-parity violation at the two-loop level jJl]. We are going 
to concentrate here on models with vectorlike tops partners. In these models, the effects on 
the Higgs mass have been so far just studied in the effective potential approach at one-loop. 
Also a careful analysis of the threshold corrections to the standard model (SM) gauge and 
Yukawa couplings has been not performed to our knowledge so far. However, it is well known 
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from the MSSM that the SUSY threshold corrections and one-loop momentnm dependent 
effects can alter the Higgs mass by several GeV j33]. Of conrse, also two-loop corrections 
involving colonred states are crncial in the MSSM and it wonldn’t be possible to reach a 
mass of 125 GeV withont them [56H69]. As soon as the Yukawa-like interactions of the new 
(s)tops become large, one should expect that effects of a similar size than in the MSSM 
sector appear. Therefore, we make a careful analysis of all three effects: we calculate the 
full one-loop threshold corrections to get an accurate prediction of the running gauge and 
Yukawa couplings at the SUSY scale, we include the entire dependence of external momenta 
at the one-loop level, and we add the all two-loop corrections which are independent of 
electroweak gauge couplings. In this context, all calculations are performed within the 
SARAH [701175] - SPheno [76l [77] framework which allows for two-loop calculations in SUSY 
models beyond the MSSM [78l [79]. The obtained precision is comparable to the standard 
calculations usually employed for the MSSM based on the results of Refs. [65] - 169j . 


Finally, we extend the particle content to have a complete 10 and 10 of SU{5) in addi¬ 
tion to the MSSM particle content to get a model which is consistent with gauge coupling 
unihcation. This model has already been studied to some extent after embedding it in min¬ 
imal supergravity or GMSB [^ [5U1 [ST]. We choose here the variant where SUSY breaking 
is transmitted via gauge mediation and check for the hrst time for the hne-tuning in re¬ 
gions which are consistent with the Higgs measurements. We show that this gives usually 
a hne-tuning which can easily compete with other attempts to resurrect natural GMSB by 
including non-gauge interactions between the messenger particles and MSSM states [821 - 192] . 


This manuscript is organized as follows. We hrst introduce the minimal SUSY model with 
vectorlike top partners as well as the UV complete variant embedded in GMSB in sec. 


In sec. HH we summary briehy the main features of the tree-level masses before we explain 
in large detail the calculation of the one- and two-loop corrections. The numerical results 
are given in secs. [IV[ and |V} In sec. 


IV we discuss the impact of the diherent corrections at 


one- and two-loop on the SM-like Higgs mass using a SUSY scale input, before we analyse 
in sec. |V]the hne-tuning of the GMSB embedding. We conclude in sec. 


VI 
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II. THE MSSM WITH VECTORLIKE TOPS 


A. The minimal model 
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We extend the particle content of the MSSM by a pair of right-handed vectorlike quark 
superfields T' and T'. The particle content of the model and the naming conventions for all 
chiral superfields and their spin-0 as well as | components are summarized in Tab. 


IIA 


In 


addition, we have the usual vector superfields 13, W, G which carry the gauge bosons for 
U{1)y X SU{2)l X SU{3)c as well as the gauginos \b, \w, ^g- The full superpotential for 
the model reads: 


W = UEjHd + 






( 1 ) 


Here, we skipped colour and isospin indices. The Yukawa couplings Y^, Y^ and Y^ are in 
general complex 3x3 matrices. The new interaction Yt' is a vector, but we concentrate only 
on cases where the third component Y^ has non-vanishing values. To simplify the notation, 
we define therefore 

y,? = y, (2) 

When we speak about the top-Yukawa coupling Yt, we refer to Y^^. 

The dimensionful parameters in the superpotential are the /i-parameter known from the 
MSSM, as well as the mass term M^' for the vectorlike top quark superfields, and a bilinear 
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term mf mixing the new states and the MSSM ones even before electroweak symmetry 
breaking (EWSB). 

The soft-SUSY breaking terms for the model are 

-£ ={Tfke,H^ + + Tiq.i'H^ + BtI'I' + Blu,F + h.c.) 

+ + m\ijd*dj + mlijq*qj + + mU*lj + 

+ ?7ip|t I + I T h.c.') + (iffiAijAs + M 2 \w^w T M^Xg^g T h.c.) 

(3) 

In general, the T- and B- parameters are complex tensors of appropriate dimension, while 
the mass soft-terms for scalars are hermitian matrices, or vectors or scalars. The gaugino 
mass terms are complex scalar. However, we are going to neglect CP violation in the soft- 
sector, i.e. all parameters are taken to be real. For the trilinear soft-term of Yf we use a 
similar short-hand notation Tj) = Tt' in the following. 


B. UV completion and fine-tuning 

1. Gauge coupling unification 

If we just include the right-handed top superhelds, the model is not consistent with gauge 
coupling unihcation. To cure this problem, additional helds have to be added. The minimal 
choice is to add a pair of complete 10-plets under SU{5) which contain the states we are 
interested in, but also vectorlike left-handed quarks {Q', Q') and vector-like right-handed 
leptons {E', E'). To generate mass terms for all components of the 10 and 10, the following 
extension of the superpotential is needed: 

AlP = Mq,Q'Q' + Me'E'E'. (4) 

Here, the Q-helds have quantum numbers (^,2,3), (—1,2,3), while the vector-like leptons 
E\ E' carry quantum numbers (±1,1,1) with respect to U{1)y x SU{2)l x SU{3)c. We 
are going to assume that no further interactions between these additional states and the 
MSSM sector are present, i.e. these particles are only spectators when calculating the SUSY 
mass corrections. Nevertheless, because of their impact on the SUSY RGBs and also on the 
threshold corrections to the SM gauge couplings they can play an important role. We can 
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see this already at the one-loop RGBs of the gauge couplings for the minimal model and the 


UV complete version: 

= (y + 

(5) 

= (1 + ^^uv) gl 

(6) 

= (“2 -k 25uv) gh 

(7) 

where we parametrized the /3 function as 

B — ^ /^(i) + ^ fl(2) 1 

^9^ - lQ^,2^9^ (167r2)2^3» ^ 

(8) 


For 6uv = 0 we obtain the minimal model, while 6uv = 1 describes the UV complete version. 
In Fig. [^the re-established gauge unification can be observed. The one-loop f3 functions of 



FIG. 1. Running of the gauge couplings z = 1, 2, 3, at 1-loop). The dashed lines belong to 

the minimal vectorlike top model and the full lines to the UV-completed model. The dotted lines 
represent the SM-only running up to MsusY = 1500 GeV. 

the Yukawa couplings are the same in both model variants and read 

ki = n(3n'n + Fn + 3Tr(y,y,t) - - 3gi - + Tr(y;n')) + (nr;) 

(9) 

df = 3nn'n + n(3Tr(yjyj) - 3g| - gg? + Tr(ny;)) (10) 

Y'.. = + 3Tl (3K„Fj) + YJy; + 6(y,,y;) - llgj - 3gl - “g^)n,i, (11) 

Y’ = 3y,.,„(nr;)+ r„(3y;y„ + rjy, + 3iv(r.y;) + 3(r,-r,;) - (Ig? - sgl - “gf) 

( 12 ) 
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Yf 


FIG. 2. This plot shows the scale Me at which the Landau pole arises as a function of 1^/. The red 
lines are for the minimal model, the blue lines for the UV complete version. For the dotted lines 
we used tan j3 = 10, for the full ones tan (3 = 60. 


We can use these RGBs to make a quick check for the cut off-scale of the theory in 
the limit of very large Yf. For this purpose, we fix at Msusy = 1-5 TeV the SM gauge 
couplings as Qi = (0.47,0.64,1.05), and consider only third generation Yukawa couplings 
Y^'^ = \/2/2AQ ■ (1.8/ cos /3, 2.4/ cos /?, 160/ sin (3) with j = e, d, u. Of course, this is a very 
simplistic setup missing many details like two-loop effects in the running and threshold 
corrections. These effects will be included in our numerical analysis. Nevertheless, one can 
already see in Fig. |^that the cut-off scale Me at which the Landau pole arises, given as a 
function of Yf, is pushed towards higher scales in the UV complete version. 


The additional soft-terms which appear because of the extended particle content are the 
following; 


— AT = rug/ 1 


|2 I 2 I 


12 I 2 I -/12 I 2 I 

I + I + 


-/\2 


+ + mYq'q' + h.c.) 


(13) 


We can now embed the UV complete version in a constrained setup to relate the SUSY 
breaking parameters. We are going to choose the setup of gauge mediated SUSY breaking 
(GMSB) which we introduce now briefly. 
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2. Gauge mediated SUSY breaking and boundary conditions 


The mediation of the SUSY breaking from the secluded to the visible sector happens in 
GMSB by messenger particles charged under SM gauge groups. The minimal model provides 
a pair of 5-plets under SU{5) which don’t have any interaction with the MSSM sector but 
due to the gauge couplings. The necessary ingredients to break SUSY are the interaction of 
the messengers, called <h, and a spurion field S described by 

W = . (14) 


S' is a gauge singlet and acquires a vacuum expectation value (VEV) along its scalar and 
auxiliary component due to hidden sector interactions, which we leave here unspecified 


(S) = M + , 


(15) 


The coupling A of eq. (14) can be absorbed into the redehnitions of M = AM and F = \F. 
With these conventions, we hnd that the fermionic components of the messengers have a 
mass M, while the scalars get masses 


0+,- = ^ ± (j)M^ , rn+_ = \/M‘^±F . (16) 

This gives the condition M^ > F. The soft breaking masses of the MSSM helds are gener¬ 
ated via loop diagrams involving the messenger particles. The gauginos receive masses 
at one-loop level while the scalar masses are generated at the two-loop. The leading 
approximations for the soft breaking masses are 

= (17) 

r=l 

ai{t) = g‘f/{4:Ti) are the running coupling constants at the scale t and Cr is the Casimir of 
the representation r. The SUSY soft breaking scales and As depend on F and M as 
follows: 

with 


9 4 


,, , , lla:^ 319x® 

/(x) ~ 1 -I- 

’ 36 450 11760 


-I- 0(x®) . (19) 








It is convenient to define 


A = 


F 

M 


( 20 ) 


For F this leads to Ac = A^ = A. Applying the general results to our (UV complete) 

model, we have the following boundary conditions at the messenger scale M for the scalar 
soft masses 




+ ^^^ 2 ) A; 

Kn = + ^92 + A 

2 2 f ^ A , ^ A \ X ‘ 
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2 2 ” 4 A 2 

^e,jj = ^e> = =-^9lAs 


’"ij; = a; 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


with j = 1,2,3. All off-diagonal entries are staying zero at the messenger scale. For the 
gaugino mass terms, we have the MSSM results 


Afi — gkc 


while all other soft-terms vanish up to two-loop 


T. =0 

Bx =0 


= Bt =0 


X = d,u, e, t' 

X = Q\E'X 


(26) 


(27) 

(28) 

(29) 

(30) 


Furthermore, we assume that the bilinear mass terms for the vector states unify at the 
messenger scale 

Mj'/ = Mqi = = Myi (31) 

We make no attempt to explain the size of /i or in this setup. There are several proposals 
how these parameters receive numerical values needed for phenomenological reasons [931 - I95] . 
We take it as given that one of these ideas is working and calculate the /i and B^ from the 
vacuum conditions. Similarly, we are also agnostic concerning the cosmological gravitino 
problem usually introduced in GMSB by the Gravitino LSP and possible solutions for it 
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[96HTn2]. 


Thus, our full set of input parameters in this setup is 

M, A, tan j3 , My , Yf 


( 32 ) 


3. Fine-tuning 


Fine-tuning addresses the question how to quantify if a model and a particular parameter 
point is natural or not. For this purpose, different measures are proposed to calculate the 
fine-tuning (FT). We are using the measure for the electroweak fine-tuning introduced in 
Refs. |1U3111U4| 


Apt = max Abs [Aq,] , 


_ 5 In M| a (9M| 
“ 9 In a M| da 


(33) 


In this setup, the sensitivity of the Z mass on the fundamental parameters at the UV scale 
is calculated, a is a set of independent parameters at this scale and gives an estimate 
of the accuracy to which the parameter a must be tuned to get the correct electroweak 
breaking scale HE]. The smaller Apt, the more natural is the model under consideration. 
We use the messenger scale M in GMSB as a reference scale and calculate the FT with 
respect to 

a = {A, My, Yp, Yt, 5 - 3 , /r, (34) 


The practical calculation of the FT in our numerical calculation works as follows: we vary 
these parameters at the messenger scale M and run the two-loop RGBs down to the SUSY 
scale. At the SUSY scale, the electroweak VEVs are calculated numerically using the mini¬ 
mization conditions of the potential and the resulting variation in the Z mass is derived. 


III. THE MASS SPECTRUM OF THE MINIMAL MODEL 

To get a good estimate of the fine-tuning by including the Higgs constraint, it is necessary 
to reduce the theoretical uncertainty of the Higgs mass prediction. Our aim is to get the 
same uncertainty as for the MSSM, namely to consider the Higgs mass in the range 

mh = (125 ± 3) GeV (35) 
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This precision can only be reached if a full one-loop calculation is done, and the dominant 
two-loop corrections are included. Since this has not been done before in literature for the 
considered model, we discuss our calculation of the mass spectrum, in particular of the 
threshold corrections and two-loop Higgs corrections, in detail. 


A. Tree-level properties 


When electroweak symmetry gets broken, the neutral Higgs states receive VEVs Vd and 
Vu and split in their CP even and odd components: 


Hd —>■ {4>d + io'd + Vd) , —)■ ((j)u + icTu + Vu) ■ 


(36) 


We have tan /d = ^ and v = \/v‘^ + v"^ ~ 246 GeV. Using these conventions, the tree-level 
mass matrix squared for the scalar Higgs particles is the same as in the MSSM. It reads in 
the basis (0^, 0^) 

I [af + ^ 2 ) - vl^ + + |/i|2 -1^91 + atjvdVu - 

(f/? + gfjvdVu - -|(^? + 9^) ( - + vj^ + + 


2,(T) 

= 



This matrix is diagonalized by Z 


H. 




(38) 


Two of the parameters in this matrix can be eliminated by the tadpole conditions for EWSB; 
dV 1 1 


Td = 


50. 


d(pd 

dV 1 


= --Vu(^Bf, ++ :^[9l +9l)vd(^-vl + vij + + |/^p) = 0 (39) 

= + afjvui^ -'Vd + vl^ - Vd'^I^B^^ + Vu(mj^^ + I/ip) = 0 (40) 


We are going to solve these equations for the squared soft-masses and when we 
consider a SUSY scale input. That leaves three free parameters in the Higgs sector at tree- 
level: tan0, /i and The last one is related to the tree-level mass squared Mp of the 
physical pseudo-scalar via 

(41) 


Q _ _ 

^ tan 0-1-1/ tan 0 


Ml 


However, when we consider the UV completion, m‘jj and m‘jj are fixed at the SUSY scale 


and we are going to solve the above equations (39) and (40) for /x and B^. 


11 








Also, the mass matrices for the CP-odd and charged Higgs bosons, for down (s)quarks, 
charged and neutral (s)leptons, as well as for neutralino and charginos are identical to the 
MSSM. Only in the up (s)quark sector things change because of the additional top-like 
states. The scalar mass matrix that links the left- and right-handed MSSM up-squarks and 
the new vector-like states is given in the basis of by 


72 I (2 (^MT'm*, + mT j + vlY*■ 

with the diagonal entries 


(42) 


mu^ui = - ^ ( - 3^2 + ^?) 1 ( - ^ ( 2 mJ + vl (Y;Yt> + FJk) ) (43) 

rriu^u*^ = ^(2(m;,mf + j + 1^2^^_ ^2 ^ (44) 

^i'Y = \ - vl + vfj (45) 

= (ml + iMrf + +^9l(^-vl + (46) 


This matrix is diagonalized by Z^: 




m 


dia 

2^u 


(47) 


and we have eight mass eigenstates called Ui in the following. Similarly, in the fermionic 
counterpart we choose the basis {uL,i,i'*) / {u*R,id'*i 3 ^ ■ The mass matrix in this basis reads 


wsVnHj -^VuYf 


m„= I 


mt' 


M 


(48) 


T' 


Here, we need two rotation matrices and to diagonalize this matrix. 


Ur m, 


.U 7 ' = m 


dia 


(49) 


The four generations of mass eigenstates are called Ui where the hrst three generations 
correspond to the up, charm and top quark. 
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B. Calculation of the Higgs masses at one- and two-loop 


In this section we give details about the calculation of the Higgs masses at the one- and 
two-loop level. We have performed all calculations with the combination of the software 
packages SARAH and SPheno which automatize all relevant steps. There are three changes 
compared to the calculation of the Higgs masses in the MSSM: 


1. The new vectorlike states change the threshold corrections at Mz to derive the gauge 
and Yukawa couplings in DR scheme from the measured SM couplings and fermion 
masses. SARAH and SPheno applies and generalizes the procedure of Ref. |55] to make 
this matching. We give more details about the main differences compared to the 
MSSM in sec. IHIBll 


2. At the one-loop level new contributions of 0{ati) arise. These corrections are widely 
discussed in literature and are known to be able to give a push of many GeV to the 
Higgs mass. While these corrections so far have just been calculated in the effective 
potential approach, SARAH and SPheno perform the full one-loop corrections in a dia¬ 
grammatic way including the dependence of the external momenta. This calculation 
is again a generalization of the renormalization procedure presented in Ref. [55]. We 


explain this calculation and the difference to the MSSM more detailed in sec. HIB 2 


3. At the two-loop level, new corrections 0{at'{as + at + at + at')) arise. The impor¬ 
tance of these corrections was unknown up to now. However, with the generic results 
of Ref. |106] for the two-loop effective potential implemented into SARAH |78| . a nu¬ 
merical derivation in analogy to Ref. HE] allows to obtain the two-loop self-energies 
at vanishing external momentum for the scalars which get a VEV. Moreover, since 
Ref. [^, a fully equivalent and diagrammatic calculation in the limit = 0 can also 
be performed by SARAH and SPheno. Both approaches are used to cross-check the 


two-loop results. We give more details about this calculation in sec. HI B 3 


1. Threshold corrections 

The presence of additional vectorlike states change the relations between the running DR 
parameters and the measured SM parameters. In the gauge sector, the relation between the 


13 










SM couplings (MS scheme with hve flavours) and the DR ones are 


a^\Mz) 


1 - Aa{Mz) ’ 

aP’^^{Mz) 

1 — Aas{Mz) 


(50) 

(51) 


Here, and are taken as input and receive corrections from the top loops as 

well as from new physics. For the minimal model, the thresholds read 


=£ (i - ¥ Eta log - i Eti log + Aa“==^“(A)) (52) 

AaM = i (-1 Eta log - I Et. log ^ (53) 


We absorbed all corrections which don’t change with respect to the MSSM in AQ;|^®®^(/i) 
and AQ;’^®®^(/i). Note, this does not include the up-squark sector, now consisting of 8 
squarks, to prevent double counting. In the case of the UV complete model, additional 
terms of the same form show up. 

To relate a to the running couplings gi and g 2 , the running Weinberg angle sin© and the 
electroweak VEV in DR scheme are needed. Also here the vector-like tops enter because 
of the new loop corrections to the mass shifts (5M| and SM^ of the gauge bosons. The 
corrections from the extended (s)top sector to the transversal self-energies are 
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(55) 


14 



with 


(56) 

(57) 


Ho{p,mi,m2) =4:B22{p,mi,m2) + Go{p,mi,m2) , 

Go{p,mi,m2) ={p^ - ml- m^)So(p, mi, m2) - Ao("ii) - ^o("i2) , 



(58) 


The appearing vertices are given in appendix A 1 All other contributions are identical 


to the MSSM and given for instance in Ref. [HS]- With that information, v and sin^ are 
calculated by 



(59) 


(60) 


Here, is the Fermi constant and 6r doesn’t receive new corrections compared to the 
MSSM (Expressions for 5^ can be found in |108| ). Also here the spectator fields in the UV 
complete version will show up in a similar way because their contributions don’t vanish even 
in the limit that all superpotential and soft-breaking interactions of those are assumed to 
vanish. 

The running Yukawa couplings are also calculated in an iterative way. We concentrate on 
the quark sector, because the leptons don’t get new contributions from the new vector-like 
quarks at one-loop. This is also true for the UV complete model because these contributions 
are proportional to the superpotential interactions which we assume to vanish for the E' 
and Q' states. The starting point are the running fermion masses in DR obtained from the 
pole masses given as input; 
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The two-loop parts are taken from Ref. |in9[ 1110] , The DR masses are matched to 
eigenvalues of the loop-corrected fermion mass matrices calculated as 


= mf^ - - T.R{pi)m)"’ - ’^Lipf) 


(T), 


(65) 

( 66 ) 

the 

(67) 


Here, the pure QCD and QED corrections are dropped in the self-energies S because they 
are already absorbed in the running DR masses. The self-energy contributions from the 
extended (s)top sector to down-quarks are 

2 4 


TTl T ^ 


a=l 6=1 
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( 68 ) 
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6=1 
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L* j-^L 

dj,W~,Uf, di,W~,Uf, 
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(69) 

(70) 


The full self-energies in the up-quark sector read now 


m rf? 

Ui,"'Ub-‘- Ui,ha,Uf, 


Kf (p') = ^0 (p^ m^^, m^_) + Bq [p^, J 

+ m--Ro (p^ ml ., m|)rf%_ ,-^rJ + m.^Ro (/, <, niAo)r^;,„,,^or2 
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- - 4:Bo(^p\ml^,ml^T^*^z,u,rnu,T^^,z,u, ( 71 ) 
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Kt(p") = Kfip") 


\U,R/ 2' 


(LoR) 


(72) 

(73) 


Because of the length of the expressions eqs. (71 73), the sums over internal generation 
indices a and b are understood. All necessary vertices are listed in Appendix |A2|f[ The 
eigenvalues of rrLj^\p‘f) must fulfill 

n^o- — .rr,2 'll _ ^DR,SM ^DR,SMn 

Big m^ {p -m^J - [^d ) 


Eig K«(y = ml)] = (m™.™, ™f’'-“.-»r) 


(74) 
(76) 

with the DR-masses taken from eqs. (61 63). In addition, the rotation matrices diagonaliz¬ 
ing and rriu^'^ are constrained by the measurement of the CKM matrix. One can use 

these conditions and invert eq. (67) to get expressions for the tree-level mass matrices, which 
are then used to calculated and Since the self-energies depend on the Yukawa 

matrices, the entire calculation has to be numerically iterated until a stable point is reached. 


After the calculation of the gauge and Yukawa couplings at Mz is finished, the two- 
loop RGBs shown in Appendix are used to run the couplings up to Msusy- Since in all 
calculations the masses of the SUSY states at Mz are needed, also a two-loop running of all 
parameters from Msusy to Mz is done to get the running tree-level masses at Mz- 

The effect of the threshold corrections on the running value of the top Yukawa coupling 
(Yp^) at the SUSY scale as a function of 17/ is shown in Fig. We have used two different 
^ the rotation matrices of the external states (marked as i in the expressions for S) have to replaced by 
the identity matrix since the corrections to the mass matrices are calculated 
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Yt> 


FIG. 3. Running top Yukawa coupling {Y^^) at the SUSY scale as function of Y^/ for two different 
values of 1.0 TeV (blue) and 3.0 TeV (dotted red). 

values of Mt''- 1 and 3 TeV. In addition, we fixed tan/3 = 3 and all soft-masses to 1.5 TeV. 
In total, this effect can be as large as a few percent and is larger for smaller Mt' because 
the t — t' mixing becomes larger. This already gives an important change in the MSSM-like 
corrections to the Higgs states which turn out to be of order of a few GeV, as we will see. 
One might wonder why the values for the top Yukawa don’t agree for T)/ = 0. The reason 
is that the threshold corrections to are always present and they depend on M^/, even if 
other couplings of the vectorlike states are absent. This changes the prediction for which 
enters (i) the SM and MSSM part of the thresholds corrections, and (ii) the RGBs when 
running from Mz to Msusy- 



Yt> 


FIG. 4. Absolute size iV'gl of the GKM entries between the vectorlike top states and the SM down 
quarks q = u,s, b. The colour code is | V)'d| (full blue), |Vt's| (dotted red), and |Vt'b| (dashed green). 
We fixed here M^' = 1 TeV. 

While a study of flavour physics in this model is beyond the scope of this paper, we 
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want to briefly comment on the expected effects. The CKM matrix in this model is a 4 x 3 
matrix and we adjnst the Ynkawa conplings and in onr stndy in a way that the 3x3 
snb-matrix assigning the conplings between SM-qnarks is in agreement with measnrements. 
The last colnmn of the CKM matrix carries the elements Vfq which define the size of the 
flavonr changing charged cnrrents between the vectorlike top and the SM down-qnarks. The 
size of \Vt'q\ is constrained by the measnrements of flavonr violating processes which are 
known to a high precision and which are in agreement with SM predictions. In Ref. m 
the following limits were derived at 3 ct: 

\Vt,d\ < 0.01, < 0.01, \Vt^b\ < 0.27 (76) 

We show the prediction of these elements as a fnnction of T)/ in Fig. for Mt' = 1 TeV. 
One can see that the obtained valnes are well below the cnrrent bonnds. The main reason 
for this is that we assnme Y)/ and Y)? to vanish. 


2. One-loop corrections 


A generic one-loop calcnlation with SARAH and SPheno was introdnced in Ref. ma. The 
procednre for this is as follows. First, all rnnning tree-level parameters are calcnlated at 
the SUSY scale. The gi {i = 1, 2, 3) and U {i = e, d,u) are obtained by rnnning np the DR 
valnes calcnlated at Mz, the Higgs soft-masses and are derived from the tadpole 


eqnations eqs. (39)-(40). Using these valnes all tree-level masses are obtained and 5M| is 
calcnlated. This qnantity is needed to get the correct electroweak VEVs at the SUSY scale 
from the Z-boson pole mass and tan /3 via 


,SUSY 




cos/3, = n®^®^sin/3 (77) 


With these valnes the tree-level masses are re-calcnlated and the calculation of the one-loop 
corrections is started. Here, first the one-loop corrections to the tadpole equations Ti 
are needed. The changes compared to the MSSM stemming from vectorlike tops are: 

4 8 


yi,u^Y) _ 


+ 65/A 




m,. 




+ F 
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i^Ua j'Rq, 




0 ’"s 


(78) 


a=l 


a=l 


with i = u,d. All other corrections are identical to the results of Ref. |55| . Afterwards, 
we need the one-loop corrections to the scalar Higgs mass matrix. Here, the vectorlike top 
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quarks contribute to the scalar self-energy 11 (p^) 




a=l b=l 


3J2^0 (< + 3 ^ ^ So (p\ml^,ml) J, 


a=l 


a=l 6=1 


Pj j'^a 


(79) 


The necessary vertices to calculate and 11*’*'(p^) are given in Appendix A3 We can 

now express the one-loop corrected mass matrix of the scalar Higgs by 

m^’*'^'^^(p^) -I-n“’'“(p^)-I- j 


n 


MSSM,'_2 


{P^) + 


Vu 

J^^MSSM^(l) 
h’d 'd,fi ^ 


0 

J_^MSSM^(1) 

Vu ^ 


(80) 


Here, Hg^M are the MSSM results without any contributions from up 

(s)quarks. The eigenvalues of correspond to the loop corrected Higgs masses. 

Since, m^’*'^'^^(p^) is a function of the external momentum, this calculation is usually iterated 
until a stable solution for each eigenvalue is found. 

Previously, the one-loop corrections in this model have been calculated in the effective 
potential approach |46) . This calculation is equivalent to ours in the limit p^ —)■ 0. Thus, 
by checking this limit we can easily estimate the error introduced in these calculations by 
that approximation. Since the additional fermions and the scalars are usually heavier than 
the desired Higgs mass of 125 GeV, one can expect that the momentum effects are rather 
moderate. However, before we discuss this in detail, we go one step further to the two-loop 
corrections. 


3. Two-loop corrections 

It is very well known that two-loop corrections in the MSSM are crucial: they can give a 
large push to the Higgs mass and are the only chance to get agreement between the Higgs 
mass in the MSSM for moderate SUSY masses (< 2 TeV) and the measurement of about 
125 GeV. This mass is out of reach only using one-loop corrections. This is not necessarily 
the case for models with vectorlike quarks: if the new couplings to the SM-like Higgs are 
large enough, even one-loop corrections might be sufficient to hnd a sufficiently large Higgs 
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mass. Nevertheless, there are good reasons to consider also the two-loop corrections: to 
be able to make any meaningful statement in the considered model if a point is excluded, 
the difference to the measurement must be larger than the theoretical uncertainty. At one- 
loop the theoretical uncertainty in the Higgs mass prediction can easily be 10 GeV or more, 
i.e. it is not possible at all to restrict many regions of the parameter space by a one-loop 
calculation. Of course, also the opposite might happen; points which are in good agreement 
at one-loop can be ruled out by a two-loop calculation. 

For this reason, we are going to give details about a two-loop calculation including the 
dominant corrections. ’Dominant’ in this context means all contributions excluding those 
of the electroweak gauge couplings gi and g 2 - That’s the same precision which is also 
usually considered for the MSSM. The remaining electroweak corrections, together with the 
missing momentum dependence and the unknown higher-order corrections are estimated to a 
remaining uncertainty of about 3 GeV. In the MSSM the most dominant two-loop corrections 



FIG. 5. Two-loop diagrams giving contributions to the effective potential 0{atas) and 0{atias)- 
Here, the indices of up-quark generations (ttj) run from 1 to 4, and those of up-squark generations 
{ui) from 1 to 8. 

are those involving the strong coupling constant g^ because of large colour factors. The 
diagrams which contribute in the MSSM are depicted in Fig. In the model at hand with 
vectorlike tops, the diagrams are actually the same but with a sum over a larger number 
of (s)quark generations. The obtained corrections from these diagrams are 0{atas) and 
0{at'as) with at = = {Y^^f/ATr. 

The next important contributions from the MSSM are those of 0 (q:^). These come from 
diagrams involving (s)tops and Higgs states respectively Higgsinos. Also here, the diagrams 
shown in Fig. are the same as in the MSSM, but the sums over (s)fermion generations 
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FIG. 6. Two-loop diagrams giving contributions to the effective potential 0{al), 0{a^,), and 
0{atat>). Here, = {h,H,G^,A^}, <!>=*=}. The index ranges are: 

ch(l, 2); x°(l - 4); 2); u{l - 4); d(l - 3); h(l - 8); d(l - 6). 


are extended. These diagrams give contributions of the order 0{a^), O^afat) and 0{a^,). 
Also the corrections 0{at{ab + ar)) with at = (Y'^^^)^/(47r), ar = (T,-^^)^/(47r) are known 
in the MSSM. Especially for moderate values of tan f3 these corrections are less impor¬ 
tant. Nevertheless, in our calculations also these corrections together with the counterparts 
0{at>{aii + Ur)) are included. 

SARAH and SPheno offer two possibilities to calculate the two-loop corrections to scalar 
Higgs masses: either a purely effective potential calculation can be done. In that case, the 
diagrams as shown in Figs. and are calculated to get and the derivatives of 

the results with respect to the Higgs VEVs are taken to get the two-loop corrections to the 
tadpoles and self-energies 


, (2L) Q2v^U2L) 

' dvi dvidvj 


(81) 


However, this involves a numerical derivation which sometimes suffers from numerical prob¬ 
lems and rather large uncertainties. Thus, the second method implemented in SARAH and 
SPheno is often the preferred one: this method employs a diagrammatic calculation where 
the external Higgs legs explicitly show up. Even if this leads to a much bigger set of two-loop 
diagrams, the calculation is not necessarily slower. All diagrams are evaluated in the limit 
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—)■ 0, i.e. the results give equivalent results for and as the first method does. 


Given the two-loop corrections, the loop-corrected Higgs mass can be expressed by 


(p^) = (p^) -|- (0) 






(2L)^ 


— (6t 

Vi, ^ 


(IL) 

u 




( 82 ) 


Here, we have no longer distinguished between corrections involving vectorlike tops or not, 
but used H^^^) and for the sum of all contributions. The eigenvalues fulhlling 

Eig(m^’^^^^(m^J) = m\, are associated with the scalar pole masses. In the following, the 
smaller value corresponds to the SM-like Higgs boson and we are going to use the short 
notation rrih = for h- 

Before we turn to the full calculation, we want to discuss briefly the importance of the 
different contributions at two-loop. For this purpose we depict in Fig.j^the different two-loop 
contributions to the Higgs mass matrix: 


Uij = t = d,u (83) 

It turns out that the corrections 0{at'ab) are negligible. The corrections 0{at/ar) are even 
much smaller and therefore not shown in Fig. We consider here two different cases: 
vanishing Tf and Tf = 2000 GeV ■ T)/. In both cases we hnd that the most dominant 
contributions are those involving the strong interaction what’s similar to the MSSM. The 
next important ones are those 0{at'at), while the 0{a^,) contributions are moderately small. 
Here, the difference compared to the MSSM corrections 0{asat) and 0{af) which often 
cancel to some extent, is that here the contributions come with the same sign. We also 
see that for most contributions the impact on the (1,1) element is the largest one, i.e. the 
dominant part of these contributions come from F-terms ~ pYf. Thus, the new two-loop 
corrections are expected to be more important for parameter regions where the light Higgs 
has a larger fraction. The main differences between the cases of vanishing and non¬ 
vanishing Tf is that the corrections involving the strong interaction to (1,1) become smaller, 
while those to the (2,2) increase. Also the 0{atiat) contributions to the (2,2) are enhanced. 
Thus, another region where the new two-loop corrections are expected to become important 
are those with large trilinear soft-terms T^/. 
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FIG. 7. Two-loop contributions the Higgs mass matrix involving vectorlike (s)quark. We used here 
Mj-/ = 1.0 TeV and put all soft-mass terms to 1.5 TeV. On the left, we set T^/ = 0, on the right 
Tf/ = 2.0 TeV • Yt'. Dashed lines are for the (1,1) element, full lines for the (2,2) one, and dotted 
lines for the off-diagonal contribution. In the first three rows we plot the individual contributions 
0(a^/), 0{at'<yt)^ ©(a^/a^), while the last row shows the sum of all contributions. 
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IV. RESULTS - PART I: THE HIGGS MASS 


Before we turn to our main results, namely the discussion of the fine-tuning in the UV 
complete model, we want to discuss the importance of the different Higgs mass corrections 
we have included. For this reason we consider first the minimal model with the MSSM 
extended by vector-like tops only. To deal with the large number of free parameters at the 
SUSY scale when not considering an UV embedding, we make the following assumptions 
about the MSSM soft masses: 

m\ = m? = = m? = m? = 1 ■ (l.STeV)^ 

Ml = 0.5 TeV, M 2 = 1.0 TeV, M 3 = 2.0 TeV 

= Td = Te = 0 

Moreover, we fix usually the MSSM parameters 

p = 1.0 TeV, Ml = (lTeV)2 

and for the new sector we assume if not stated otherwise 

Tf =mt = Bt = Q 
Trip = m% = (1.5 TeV)^ 

t u 

In addition, the most important SM parameters were chosen as 

c^{Mz) = 0.1180 , mp(mb) = 4.2 GeV, mf'" = 173.2 GeV 

As already mentioned we employ the combination of the computer tools SPheno and SARAH 
for all numerical calculations: we have implemented the minimal model with vectorlike tops 
as well as the UV complete variant in SARAH version 4.5.3 and the model files will become 
public with the next release of SARAH. SARAH was used to generate Fortran code for SPheno. 
The obtained Fortran routines include automatically all new features from vectorlike stops 
discussed in the last sections which are necessary for the precise Higgs mass calculation. Also 
routines for the calculation of flavour observables and decays widths are generated by SARAH. 
However, we will not go into details in these aspects of this model here. We are just using 
the FlavorKit results m to double check that all points are in agreement with current 
bounds from flavour observables. This is, of course, expected as we already discussed in 
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sec. 


Ill B 1 The Fortran code written by SARAH was compiled together with SPheno version 


3.3.6. For all parameter scans in the following we have used the Mathematica package SSP 

m- 


A. The difference between one-loop effective potential, full one-loop and two-loop 


We check the importance of the corrections calculated here for the hrst time. For this 
purpose we compare in Figs. |^- IT the prediction for the Higgs mass calculated 


(i) at one-loop with vanishing external momenta but including thresholds, 

(ii) at one-loop with full momentum dependence but neglecting the threshold corrections 
to SM gauge and Yukawa couplings, 

(iii) at full one-loop including the full momentum dependence and all threshold corrections, 

(iv) at full one-loop with dominant two-loop corrections. 


The one-loop calculation without external momenta is equivalent to the calculation per¬ 
formed in the effective potential approach. For all three Figures we have put M^' = 1 TeV. 
In Fig. I^we compare the results for two different values of tan/?: 2 and 10. While there is 
a large difference already at tree-level, the impact of the loop corrections is similar for both 
values of tan/?. Thus, we hnd that ruh — 125 GeV is found for W ~ 0.9 (0.6) for tan/? = 2 
(10). Including the momentum dependence in the one-loop calculation of the vectorlike 
states can account for changes up to 2 GeV for large Y^/ and are negative. In contrast, for 
the considered scenario the two-loop corrections are of a similar size, but positive. However, 
the biggest difference are caused by the threshold corrections. Since these can have a large 
impact on the top Yukawa couplings, we hnd that the prediction of the SM-like Higgs mass 
can deviate by up to 5 GeV. This effect is more pronounced for smaller tan/?. Note, even 
in the limit V/ —>■ 0, we hnd a shift by about 1 GeV compared to the calculation using 
only MSSM results. The reason is that the threshold corrections to don’t vanish even 
in this limit. Therefore, the running value of the top Yukawa coupling entering the loop 
calculations changes slightly, which has still a visible ehect on the Higgs mass. The absolute 
size of the one-loop corrections can grow up to 30 GeV for both values of tan /?, while the 
two-loop corrections are smaller by about a factor of 10. When we compare these numbers 
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FIG. 8. Top left: light Higgs mass as function of Yf. The red line corresponds to the effec¬ 
tive potential calculation at one-loop, orange is the one-loop corrections with external momenta 
but neglecting the new threshold correction stemming from vectorlike states, blue is the full one- 
loop calculation including the momentum dependence and all thresholds, and green includes the 
dominant two-loop corrections together with the full one-loop correction. Top right: impact of 
the threshold corrections (red), the momentum dependence at one-loop (orange) and the two-loop 
corrections (green), given as the difference Am^ = mh — m/i(lL,p^ ^ 0,all thresholds). Bottom 
left: absolute size of the one- (blue) and two-loop (green) corrections stemming from the vectorlike 
states. Note, for better readability we re-scaled the two-loop corrections by a factor of 10. Bottom 
right: relative importance of the one- (blue) and two-loop (green) corrections normalized to the 
size of the purely MSSM-like corrections. The full lines are for tan /3 = 10 and the dotted one are 
for tan/3 = 2. We used here Mp/ = 1.0 TeV, = 0. 

with the purely MSSM corrections, we see that the one-loop corrections can become as 
important as the MSSM ones, while the two-loop corrections can reach about half the size 
of the MSSM two-loop corrections. 
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FIG. 9. The plots show the same results as in Fig. when including non-vanishing Tf. We used 
Tf = 2.0 TeV • Yf, tan/3 = 5 and = —2500 GeV. The full lines are for Bt' = 0, while the 
dotted ones correspond to Bt' = (1.5 TeV)^. 


We have identified in sec. Ill B 3 two regions where the new two-loop corrections are 
expected to be even more important. The hrst region is the one with non-vanishing Tf/. 
This is studied in Fig. where we set Tf = 2000 GeV ■ Yf. In addition, we check also the 
effect of Bt'- For Bt' = 0 the differences to the results with Tf = 0 are not very large: the 
corrections from the momentum dependence and the two-loop terms are of the same size 
and come with different signs. The largest effect is again from the threshold corrections. 
However, if B^' becomes large and causes a mass splitting for the vectorlike stops, the picture 
changes. Now, the most important effect comes from the two-loop corrections which can 
become as important as the MSSM ones. For Yf values of 0(1) this can reduce the Higgs 
mass prediction by more than 10 GeV and easily over-compensate the two-loop corrections 
from the MSSM sector. 

The other region we identihed where the two-loop corrections can be important is the 
one where the SM-like Higgs has a larger down-type fraction. This happens if M\ becomes 
small. We discuss this case in Fig. I^for zero and non-zero B^' again. In particular for the 
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Yt> Yt, 


FIG. 10. The plots show the same results as in Fig. I for smaller M\ = 10^ GeV^. We put 
Tt' = Tu = 0, and tan/3 = 3. The dashed lines are for Bt' = 0, while the full ones correspond to 
Bt' = (1.5 TeV)2. 

large Bt' the two-loop contributions can clearly make the biggest effect compared to the 
incomplete calculations used so far. These are again negative and can reduce the SM-like 
Higgs mass by up to 8 GeV. Thus, while it seems that one can reach the preferred mass of 
125 GeV at one-loop with V/ < 1, with the two-loop corrections this is not possible for the 
considered combination of parameters. Although if Bt' is taken to be zero, the effect can 
still be large and the overall size of the new two-loop corrections is still in the ballpark of 
the MSSM corrections. 


B. Dependence on the vectorlike masses, stop masses, and the gaugino mass 


As a next step we want to understand the dependence of the loop corrections on the 
involved masses a bit more. We start with the dependence on the vectorlike mass parameter 


Mt' and Bt' and show in Fig. IT the Higgs mass at the one- and two-loop level. At one- 
loop we have the well-known picture that the corrections quickly decrease with increasing 
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FIG. 11. Contour lines of constant irih at one- (left) and two-loop (middle) in the plane. 

The plots on the right column show the size of the two-loop corrections involving vectorlike states. 
The plots in the hrst row are for 1^/ = 1.0 with T^/ = 0 and in the second for T)/ = 0.7 with 
Tt> = 1400 GeV . 


mass of the vectorlike states, while the dependence on Bj" is small and just shows up for 
smallish M^' of 1 TeV and smaller and large \Bt'\ > 2.0 TeV^ for h)/ = 1.0 and Tti = 0. 
This general picture does, of course, not change at two-loop but we hnd a shift by several 
GeV usually dominated by the MSSM-like corrections. The two-loop corrections from the 


vectorlike states are singled out in the right column of Fig. IT They don’t show this strong 
Mt' dependence as the one-loop corrections do, and actually slightly increase with larger 
Mt'- Also the dependence on Bt' is more pronounced at two-loop. If we go for smaller Tj' 
and turn on Tfi the one-loop corrections in total become smaller and are less dependent on 
Bt'- However, the sensitivity at two-loop and Mt' and Bt' is nearly the same, but just the 
total size of the corrections decreases. 

We have so far just concentrated on the dependence of the Higgs mass corrections on 
the new parameters absent in the MSSM. We want to dualize our discussion of the loop 
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FIG. 12. On the left: the light Higgs mass rrih as function of 1^/. Here, we used different values 
for M 3 : 1 TeV (red), 2 TeV (blue), 3 TeV (green), 4 TeV (orange). The full lines are the two-loop 
results, the dotted ones the one-loop. On the right the absolute size of the one- (blue) and two-loop 
(green) corrections involving vectorlike states. The line coding is dashed, dotted, dot-dashed, full 
for increasing M 3 . 


corrections by also briefly commenting on the impact of at least two MSSM parameters: the 
glnino mass parameter M 3 and the soft-mass for the left-handed stop, mg, 33 . We start with 


the dependence on the glnino mass shown in Fig. 1^ Here, we vary Yti and nse glnino masses 
between 1 and 4 TeV. At the one-loop level there is of conrse jnst a tiny impact on the Higgs 
mass. The small difference comes from SUSY threshold corrections. For Mt' = 1.5 TeV and 
3.0 TeV we find that with increasing M 3 the two-loop corrections O^asOif) become larger. 
Since they are negative, the prediction for ruh becomes smaller. However, for large Mt' the 
dominance of the corrections 0{a^,) is so large that this effect nearly doesn’t play any role. 


Finally, we check the impact of the soft-masses for the left-handed stops. The one- and 
two-loop corrections as fnnction of V/ and mg ,33 = 1, 2,3,4 TeV are snmmarized in Fig. 

We see that this parameter plays an important role at one-and two-loop: at one-loop, the 
corrections increase by a factor 2 when going from 1 to 4 TeV. At two-loop this effect is 
even more important and the corrections change by nearly a factor of 3. Interestingly, the 
one-loop corrections are larger for larger sqnark soft-terms, while the two-loop corrections 
increase with decreasing sqnark masses. 
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FIG. 13. On the left: the light Higgs mass m/j as function of 1^/. Here, we used different values 
for mq^ 33 : 1 TeV (red), 2 TeV (blue), 3 TeV (green), 4 TeV (orange). The full lines are the two- 
loop results, the dotted ones the one-loop. On the right the absolute size of the one- (blue) and 
two-loop (green) corrections involving vectorlike states is shown. The line coding is dashed, dotted, 
dot-dashed, full for increasing 771 ^^ 33 . 


V. RESULTS - PART II: THE FINE-TUNING IN GAUGE MEDIATED SUSY 
BREAKING 


We now turn to the consequence of the loop corrections for the fine-tuning in minimal 
GMSB. The intrinsic problem of minimal GMSB in the MSSM is that it predicts very small 
trilinear couplings. Thus, the only chance to enhance the Higgs mass via loop corrections 
is to go to very large values of A and M to get sufficiently heavy stops. When calculating 
the fine-tuning for this setup and demanding mh — 125 GeV, one finds that the fine-tuning 
A is well above 1000. Of course, in the presence of large loop corrections due to vectorlike 
states, the need of superheavy stops is relaxed and the fine-tuning is expected to improve. 


We show in Fig. 14 the fine-tuning in the (tan/?,y)/) plane for different constraints for the 
Higgs mass within the theoretical uncertainty: (i) ruh = 122 GeV, (ii) ruh = 125 GeV, (iii) 
mh = 128 GeV. For the vectorlike states, masses of 500 and 1000 GeV were used at the 
messenger scale. 


One finds that the fine-tuning quickly drops with increasing Yf because lighter SUSY 
states are sufficient to push the Higgs mass to the desired level. For very large Yf of 0(1) 
and the looser constraint of mh > 122 GeV, even a fine-tuning of about 100 seems possible. 
There is also another, very interesting observation: even for Yf = 0 the fine-tuning in 
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FIG. 14. Contours of overall fine-tuning A in the (tan/3, lt/)-plane demanding a Higgs mass ruh = 
128 GeV (top), rrih = 125 GeV (middle), and mh = 122 GeV (bottom) for the UV complete 
variant of the model. We fixed here M = 10^ GeV and Myt =0.5 TeV (left column), respectively, 
Myi = 1.0 TeV (right column). The red dashed lines indicate the gluino mass in GeV. 
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this model is not as bad as one expects it from the MSSM. The reason is that the strong 
interaction at the messenger scale is larger compared to MSSM expectations because of the 
different running. Therefore, for the same value of A, the squarks are already signihcantly 
heavier and lead to larger Higgs mass corrections. 
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FIG. 15. Minimal fine-tuning for given Higgs mass and gluino mass rrig. We fixed here M = 
lO'^ GeV and Mt' = 1 TeV and scanned over tan /3, Yt> and A. 

However, including the bounds from direct SUSY searches has a large impact: the points 
with a small hne-tuning are excluded because of the light gluino mass. That’s completely 
different to the GMSB variant of the MSSM where the Higgs mass pushes the hne-tuning of 
the model to higher values. In this model, the vanishing trilinear couplings at the messenger 
scale just play a subdominant role concerning the hne-tuning, but the gluino mass demands 
a larger SUSY scale A, which increases the hne-tuning. The situation wouldn’t change if we 
go to larger Messenger masses to increase the running because the one-loop /d-function of 
Afa vanishes in this model and the mass is actually slightly decreasing with increasing M. 
Moreover, it’s a general feature of GMSB that the gaugino masses are not very sensitive to 
the messenger scale because the leading dependence in the RGE running always drops out. 
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The running gaugino mass at the SUSY scale is related to the one at the messenger scale 
by the ratio of the corresponding gauge coupling at both scales: 


Mi{Q) = Mi{M) 


gKQ) 

qKm) 


- 9l{Q)^G 


(84) 


We show the minimal fine-tuning in the {rrig^mh) plane in Fig. 15 It is interesting that the 
fine-tuning for = 125 GeV can be smaller than for mh = 122 GeV and mh = 128 GeV 
when the gluino mass is sufficiently large. 




10 20 30 40 50 60 

tan/3 


FIG. 16. Contours of the overall fine-tuning A (left) and the mass of the lightest up-squark (right, 
full blue lines) and gluino (right, dashed red lines) in the (tan/3,demanding a Higgs mass 
nih > 122 GeV for the variant of the model without spectator fields. We fixed here M = 10^ GeV. 


For very large 1^/ where the FT becomes the best, the theory is not perturbative up to the 
GUT scale. Since there is a cut-off anyway in the theory, there is no real need to maintain 
gauge coupling unification by adding the spectator fields at the SUSY scale. Therefore, one 


might wonder what the FT of the minimal model is. This is depicted in Fig. 16 In this 
setup, the squarks are lighter for the same values of M and A because of the smaller strong 
coupling at the messenger scale. Thus, in general larger A is needed to increase the Higgs 
mass. This leads also to larger gluino masses. This is shown in Fig. where we compare 
the minimal value of A to get a Higgs mass larger than 122 GeV in the (tan/3, Y^/) plane for 
a messenger scale of again 10^ GeV, and the resulting stop and gluino masses triggered by 
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FIG. 17. Contours of constant A (black), the lightest top-squark mass (right, full blue lines) 
and gluino mass (right, dashed red lines) in the (tan/3, lt/)-plane demanding a Higgs mass m/j > 
122 GeV. All contours are given in units of TeV. On the left for the UV complete model, on the 
right for the model with only vectorlike tops. We fixed here M = 10^ GeV. 


this A. We find for the minimal model the following hne-tuning 

A ~ (230,275,320,380) 


(85) 


for m-g > (1000,1200,1400,1600) GeV and ruh > 122 GeV. 


VI. CONCLUSION 

We discussed the loop corrections to the light Higgs mass in the MSSM extended by a 
pair of vectorlike top quarks. We have improved previous calculations in literature in three 
respects: (i) we included the additional threshold corrections from the vectorlike states to 
SM gauge and Yukawa couplings, (ii) we added the full momentum dependence at the one- 
loop level, (hi) we calculated all dominant (i.e. excluding electroweak) two-loop corrections 
in the effective potential approach. It has been shown that the momentum effects can be 
sizeable and change the Higgs mass prediction by a few GeV. The effect from the threshold 
corrections turns out to be often more important. The importance of the two-loop corrections 
strongly depends on the considered parameter point. They are often a bit smaller than the 


36 








































































two-loop corrections known from the MSSM, but we also identified regions where they can 
be even larger. In these regions, the additional two-loop corrections can change the Higgs 
mass prediction by up to 10 GeV. We checked the impact of the presence of vectorlike 
states on the hne-tuning in GMSB. For this purpose, we extended the model by additional 
vectorlike quarks and leptons to have complete multiplets of SU{5). We found that the 
hne-tuning can be reduced signihcantly compared to minimal GMSB with only the MSSM 
particle content. Often, those regions with the best hne-tuning which are in agreement with 
the Higgs mass measurement are ruled out by gluino searches. Interestingly, we hnd that 
for heavy gluino masses the hne-tuning for heavier Higgs masses can be even better. In 
particular, for rrig ~ 1400 GeV, the best hne-tuning is found for a Higgs mass of roughly 
125 GeV. 
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Appendix A: Vertices 

1. Vector boson Vertices with vectorlike (s)tops 
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3 . Higgs Vertices with vectorlike (s)tops 


(A20) 


(A21) 

(A22) 

(A23) 

(A24) 

(A25) 

(A26) 

(A27) 


y.'^j^blk 


yR 

UiaUj 0hk 




^7^<5a/3 [UZ4 E + E E UZayu,ab) Z, 

a=l 6=1 a=l 


7H 

^k2 


3 3 






E V',aUhaUlg4 + E E KabUhaUhb 


a=l 


6=1 a=l 


(A28) 

(A29) 


39 



■p _ _ _ _ 

hiUjpu’^^ ^2 


- Sal + g?) 5 ; ZffzH, (,.,Z" - g„Z") 

a=l 

3 3 3 3 

+ 2Z^- (3V2g ^ y,;.„Z“.Z" - 3^2 ^ - 6g.. 

'^t',aXu,ba^k3+b^i2 

a=l a=l 6=1 a=l 

3 

- 2g-iviZSz" + 2glv„ZgZg - 6g„ 

a=l 

3 3 3 

-2[-3V2gY,Y. jZ-\-a Zg^Zg + 3V2Mt.Z';- ^ YgXaZg 

fe=l a=l a=l 

3 3 3 3 

+ 3V2^ zj* 5^ z'i,^j„ag + 6«„ ^ zJ-k,,. 5^ K.;,z‘lz| 

6=1 a=l a=l 6=1 

3 3 3 3 

+ Z Z%Jl^Zg,Zg + 3V2Z^- Y, Y y,l>.m,'.X,Z« 


6=1 a=l 


6=1 a=l 


+ E ^i3+c E E y:.^y.A,zg+3 v,yY. E n-..on,„.z£z| 


c=l 6=1 a=l 


c=l 6=1 a=l 


+ 29? E 4 +<. 2 ?? 3 +. ( 9 ^ 2 ?.? - g„Z|) + 3v^ ^ ZggT,,_^ZgZ, 

a=l a=l 

3 3 

+ 6g„ 5 ^ Z^:, 5 ^ Yl^Yg,,ZgZg, 


6=1 a=l 


sV^g-Z" ( ^ Z“,’-y,.,„Z'' + ^ Z"'- ^ n,„3Zg,^„) - 2glviZgi-Zg Z^ 


6=1 a=l 


+ 2gKzg*Z"z9 + 3V2MJ, z7/y,,.Z"Z‘ 


3V2^z5*5^m,-,„K„3Z"Z‘ 


6=1 a=l 


(A30) 


6-2 hjU]^ry1lj^^ 


V((-3g? + 9?) E Arz£(z"z"-z|z 


-4(35]zJi-y,,,„5^ ^t' ,b^fb^a^f2 + 3E^ EE Yt' ,aYuMZ!3+bZg Zj 


+3 E 2 S +0 E E yg^y-^MZg+bZg z‘ 


c=l 6=1 a=l 


+ 3E E Ar E n%.ll..« 32 ?£z"z" + g^ z'-Y^Zg.JzgZg - ZgZ' 


c=l 6=1 a=l 


40 



^2 ^7^)* 'Z-^ 'Z^ yU,=^ ryH ryH ryU , q ryU,^ \ I \7 |2 T'// T'i/ T' 

9l^k7 ^il^jl^lT 9l^k7 ^i2^j2^l7 ^ ^^k7 ^22^j2^i 


a=l 


Q \ ^ 'yU^^ \ 17* 17 -7-ff T'f/ 2 ryU,* ryH ryH ry 


b=l 


a=l 


2 ryU,* ryH ryH ryJJ 

9l^kS ^i2^j2^l8 


(A31) 


Appendix B: Renormalization Group Equations 

We give in the following the two-loop RGBs for the considered model. In general, the 
RGBs for a parameter X are dehned by 


d 


dt^ IbTT^'^^ (IbTT^)'^'^ 

Here, t = \og{Q/M), with Q the renormalization scale and M a reference scale. Bor a 
parameter x present in the MSSM we show only the difference with respect to the MSSM 
RGBs 

A/3^f = - i3(^)XSSM ^ 32 ) 

Where refers to the minimal model with vectorlike top quarks discussed here. The 
additional difference to the UV complete version of the model is given as 
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The calculation of the RGBs in SARAH is based on generic expressions given in Refs. |TT51 - 


1 . Gauge Couplings 
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A/?W = 0 (B8) 

= igl (B9) 
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2. Gaugino Mass Parameters 
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3. Trilinear Superpotential Parameters 
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4 . Bilinear Superpotential Parameters 


A,3™ = 39(y,y;) (B4 o) 

A,3f = +g( 20 g| + 9 ?) 4 i(y.yr) -99((y.y!'))' 

+ - 25(y,yjyyr) + - 75(y-yjyy;)) 

/^VVpm ^ _g_ 

t = ^ (i5Mt. (y.yr) +15 (y.y'm,.) - 8 ( 59 ? + gfj Mt' 

+ “^gUlMT^ + jaiMr^ - sAfr.((y.yr))" - 2 Mt' (y.y'yyr 

' _. \ _ / _ + _ / _. 


/p(2) _ I 2096 ^4 
~ “T 


(B41) 

(B42) 

(B43) 


225 ^' " " 9 JJ 

(y,y'm,,) [egl - 6Ti(yy„') - ^gj) 


(QglMT' - 6 Mt>Tt(yX^ -8(Yt,Y^mt>'^ - ‘^g^MT') 

^ = 2(Mr. (y^Y;)^ + ^ {5gl + gl)m,, 

+ 50 ^^ + 80glgl)mt',i 


(B44) 

(B45) 

(B46) 


zzo \ / 

- ^ ((20Mt^ (yt,y;) + 5(yt,y>i,) + + i5Tr(y,yoj + j j 

+ (-15.72 + i 5 Tr(y„yj) + i 5 (y,y;^ + gi) 

+ 5 [mt' (y^y.y;^^ + Mt> (yXYuY;)^ + (y^yiy.yjm*.) 



u^u^u^t' 

(B47) 


(B48) 


44 



= -f (B49) 

/3it = (B50) 

^ + ^Ogl + gfjMq, (B51) 

/?£^, = (lOg? ( 169 I + 9gfj + 25 ( 2569 I + 2889|g| + 297gj) + 289gj) Mq, (B52) 


5. Trilinear Soft-Breaking Parameters 
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- 12glM2(Y;^Y:Yt,)^ - 24 (y;T,) + 12^7^ 

- 23(y,Y;) (y^Y:T,)^ - 15Tr(F.Fj) + ^glirjY^Y,)^ 

- 6Tt(y,yJ) (tJf;F,)^ - 2 Tr(nF;) (tJf;F,,)^ + ^-gl(T^Y:Y,)^ + 6 ^ 7 ^ 

- i8(y*,F;) (t^y:y,)^ - i2Tt(yX) It^y:y,)^ - 2{yJy:yJy:t,)^ 


- a(yJy*tJy:y, 


a(y^y:yjy:tA - q(y^y:y^y:tA -a(y^y:tjy:y, 


- 8(y^y:t^y:y,)^ - - 2(t^y:yJy:y,)^ 

+ 450 9iTt',i + glglTt'^i + —gtTt',i + -^9iglTt',i + 8glglTt>^i 

+ ^^773T^^i + 277 ?(Fi,y;)Tt,, + 6gl(Yt,Y;'^n, + 16 ( 73 '“ 27((Ft,F;)) V*,,, 

- - 2a{y,,Y^Y^Y;^T,,, + ^(7iTr(F,Fj)Ti,, 

+ 16glTT(YX)Tt',i - l2(Yt>Y;^^T(YX)Tt',^ - 3 Tt(Y.Y^Y^ y^,, 


-dTrlYXYuYnT, 


(B63) 


A^^f3P, = 


^ ( -A(675glM2 + SOOg^M^ + 91gtM,yt',^^ + {Q75gl + SOOgt + 91(7t)T*,, 


(B64) 


6 . Bilinear Soft-Breaking Parameters 


= 3B,.(\',,y;) + 6li(Y;T,) (B65) 

YbI = +S;,( - I8(r,-y;ty„r,;) -6(r,.Fir,r;) - 9((ivy,;))' + |gj + g( 20 g= + »?) (r,,r;)) 
- ^fi(24gfM, + 75(y,rTjy;K,.) + 225(y,rT„’'y;y.) - iog^(y,;T,,) -2oogj(y,;T,,) 

+ 5(y.y,;) (2gfAf, + 40g|A/3 + 45(y,;T,,)) + 75{Y}YiY;,Tt^ + 225(y,;y„y,rr,.)) 

(B66) 

= ^ ( - 4 (75(72'M2 + 7gtM,) /x + ( 75772 " + 7gi) B,) (B67) 

^ (i 677 ?MiMt. + 80773 'M 3 M:r, + Bt> (i 5 {y^Y^) - 8 (5773' + 77 i') ) + 15 
+ 30Mr' (Y;Tt>^ + 30 (Y^m^Tt,'^^ (B68) 
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- l2glM2{Y,Y^me') - ‘^gl{Y,Y;iB,') + Qgl{Y,Y^B,') 
-2{Y,YiY,Y^B,)-2{Y,YXY^B,)-AMT>{Y;TjY:Y,^ 
-AMT^(^Y;T^Y:Y,^-^glMT>[Y;T,^+l2glMT^(^Y;T,^-16(^Y,^^ 
-4(f,WtJf;F,,) -4(F,WTjy:F*,) - ^gl(Y^m,T,) 

+ 12gl(Y^mt>Tt:') - AMt^ (Y^Y^Y^Tt^'^ - 4 Mt'(Y^Y^Y^T t^'^ 


-4(f]F,F,WT,,) -4(yJy.Fjmt.T,,) 

- 12Mt^(y;T,)Tt(yX) - 12(F,WT*,)Tr(Ky,t) 

+ (1048.?^ + MOglgl + AOOgI - 90o((Fi,y;))' - 225(F,.F,V,y;) 

- 225{y,Y^Y^Y;^ - 45(y,,y;) ( - Ibgl + 15Tr(F„Fj) + g^)^ 

- 12(F*,Fjm,,)Tr(y,tT„) 

+ (^^glM^Mr^ - UglM^MT' - s{Yt,Y^Bt>j - 32 Mt^ (Y.ITt^'^ - 

-12Mt'Tt(y^Tu^^ (B69) 

= -^{^{^OgtM, + 7gtM,)MT^ - {^Ogt + Tgf^Brj (B70) 

^ (i5(2Mr' + 2(T„y,tm,,)^ + 5^, + (yuVjB.'jJ 

+ 16 (bglM^ + glM^mt'^i - 8(^65(3 + 5(3j Bf^i'^ (B71) 

ISlfi'iMi + 405'35'3 (^Mi + M3 

+ ^glM^Mr' (y^Y;^ ^ - Ug^M^Mr' (y^Y;^ ^ - ^IBt' (y^Y;^ ^ + Gg^Br' (y^Y;^ ^ 

- 8Bt' (y^Y;^^ - 2(Yt,Y^Bt>'^ (^“^*0* “ 

- 4(Fjmt,Ti,) - 6Bt'Tt(YuY^^ “ 12MT.Tr( fJT„) 

- + 12(72Mt. - 16Mt. (f^.F;) (t.F;)^ - 4(FTFjmT) 

- 12MT.Tr(F,Fj) (t.F;)^ + ^< 72^3 (^F^F^m*.)^ - 12772M2 (f„F„W)^ 

- 12 (f;Tt) (F.FjmT)^ - 12Tr(FjT„) (f^F^W)^ - ^^ 7 ?(f.FJE,.)^ + 6^72 ^(f^FJE,.)^ 

- 6 (FtF;) - 6Tr(F.Fj) {YuY^B,)^ - ^77 ? (t^F^t)^ + 12772 ^ (t.F>t)^ 


+ 505(3 M3 ^ (ISIS'! + 50^3 


\B 


t' 


( 2 ) 

225 
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12 v,,k: 


) - 12Tr(K„yj) (T„y,;m,.)_ - 2BT.(y^YjY,Y;)^ 

yJt,y;) -2BT.(y„y„'y„y;) -4Mr.(y„yjr„y;) -4M7..(T„yJy,y;) 

/ i \ / i \ / i \ /I 

yjy„y,;). - 2 (y„yJy,y„'B,-)_ - 4(y„yiT,y„tm,,), - 2 (y„y.;y„yjB,-)_ 


4Mt. nyir^y;) -2b 




HTuYjy<i 


- (bOgt + 

(^(^4:5glM2 + SOglM^ + glM^Mqi — i^bg^ + SOgfg + gl^Bqi'^ 

= ^(-4(25 [U^ghl (Ms + M 2) + 256(73"M3 + 297(72"M2) + 289gtM, 
+ 5(71 (^16(73 (^1 + ^3) + 9(72 + ^2) ) ) Mq / + ^10(7i (^16(73 + 9(7^) 


+ 25 256(73^ + 288 ( 7^(732 + 297gt) + 289gt)BQ> 


I^bI, — — Be'^ 

pZ =-^gt{^M,ME^ - Be^) 


(B72) 

(B73) 

(B74) 


(B75) 

(B76) 

(B77) 


7. Soft-Breaking Scalar Masses 




= \j^9i ( “ 2Tr ^m^) + 2m| - 2m| - Tr 
r^m^) -hTr^TTig) -hTr^m^)) 

= J^9i ( - ml - m\ + ml + ml^j 




f^2,ii = +3Tr(^m^^ 

-h 8 m| + 8ml + Tr ^m^) ) 

1 . . 


+ 3m|^^ + 3m^^ + 6 Tr ^mg) + 8 Tr ( 


/\uv„ _ ^ ^2 1 

A (^ 2,11 - y^9i ( 
1 1 


\ / / 

^ 6 ^m| + m|) + ml + 

1 7 ^29 ^^99 


^3,1 = ( - 9(/im^^ - 45^2^^^ + 9glmjj^ + Abglmjj^ 

+ 82glml + 1605f|m| + l29mYYuYi: - 3017/mg*y^T 
-F 129ml (Yt'Y^1'j + 4fl'?Tr^m^) + 805f3Tr^m^) 


2 2 o o 2 2 

- 32c/imp 


(B78) 

(B79) 

) 

(B80) 

(B81) 

IbOglml 


if y; - 90mf^ (y,.y,f 


'sTrlm^ 
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+ ?>QglTi{mf^ - A5glTr(m]'^ + glTr(mfj +A5glTr(mfj + 80 ^ 3 X 1 

- 32^2Tr(^m2) - 160^^Tr(^m^) + 90mjjTr(YdY^^ + 30m]jTr(YeY^^ - 90m^^Tr 

- 60TT(YdYjml*^ - 30TrF]) - 60Tt(y^YJ ml*'J + 30Tr(^Femf*F;) 

+ 120Tr (F^Fjm^*) - 30Tr (F„mf FJ) ) (B82) 

A^^as,! = (5 (i6c/3 + 9 ^ 2 ) ( “ ( “ 36m| + 36m| - m| + m|,)) 


c^2,2 = ^ (3Tr(^m2) + 

A^^c^2,2 = ^ (ml, + m|) 

0-2,3 = ^ (2TT(mfj + m| + m| + j 

A^^a2,3 = 


(B83) 

(B84) 

(B85) 

(B86) 

(B87) 


A/ 3 W = 




t ',2 


( 2 (mL + + (™f iv),) 


+ if 


(B88) 



i 2 m?,,y,rp(ny;')y,.j -6m|y,r,Ti(y„y„')y,,j -6T;,iv(yJr„)y,-, 


6y,r ,Ti' I t:tJ I y,,., - ey; ,Tr 


- 6Y;,'Tr(^mXyu)y,'., - 6y;.Ti'(mJy.yj)y,,, + grfy,,,(mjy;)^ 


- 5(y.y;)y.,i(m5y;)_ - 3Tr(yyj)y.,,(mjy,r)_ + 

- io(y.y;)y^,(y>y)_ - 6Tr(yy„')y, j(yjmjf,)_ + ^gjy,r,(mfy.) 

- syy (y,y;) (mfy.), - 3yyTi(yy„t) (m 


y. -SmlY, 
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8 m],Y;jYjY:Y,, 


4 mlY;jYj^Y:Y,, 


2 mJy; W^y;y,, 


— 4( Yjm 


Y!Y*YA -4TY Y!Y*Th] -4Y4Ay:t:Th] -4TY T4Y*Y, 


- 4 y;At^t:y, 


',,) - 2 y,,,j(mJyJy„r,;), - 4 y,,j(yJm;y„yr)_ - 4 y.,,(ytym 5 yr)_ 




2 Y;AmfYAY:Y,, 


4 Y;AYAml*Y:Y,, 


4 Y;AYAY:mfY, 


+ 4 g;M;( 29 ?M,i + 5 y;,( 2 M,y.,j -y,^)) + gg?r;,r,,j - ior,;,,(y.y;)r,.g 


- loy;, yr; y^ - ey; ,Tr y.yj y,,, - ey; yr yy„' yg - 4 yjyy, y,,, 


u^uY^i j Tyj 


^ mA 


AAA 

^mf 


= ^l( 675 gt\M 2 \^ + + S 00 gt\Ms\‘ 

144 

= 

= —l( 75 gAM 2 \‘^ + 

144 


(B89) 

(B90) 

(B91) 

(B92) 


UVo{2) 


A^^ A. 


aaA 

I m"_ 


- 6((mi^ + + ml) (y,,yJy,y;) + (yrjyyr) + (y,mjyiyy:) 

+ (yyjmlyy;) + (y.yiy<,mly:) + (yy'yylmL--) + (y'yyy.) 

+ (yjyyry.) + (yi'yy.y;)) (B93) 

12 fi 

= + (B94) 

= 6((mi„ + m|) {y,Y;) + (Fymjy;) + + (t^T;)) (B95) 

= -36(m^„ +m?) {{YfY;))" + ^glM*(A0M,(Y,Y;^ + ISg^M, - 5(y;Ty)) 
+ ^ + 40glmjj^ + 2glml + 40^|m| + SO^^lMgp - 4b(YtmlYA^ 




-45 TyT, 


“ 5 ( “ + ^i) {Yt'mlY;;^ - 4gl (Yt^Y^mY^ - 80c/| (Yt>Y^m 

+ i 5 m^^(yyyiy,y;) + i 5 mi„(yyyiy,y;) + ism? (y^y^y^y;) 
+ 180m^^ (Yt/Y^YnY^'^ + 90m?, (Yt/Y^YnY^'^ + 15 (YfT^TdYA'^ 

+ 9 o(yyTiT,,y;) + i 5 (Y,mlYjYdY;^ + 9 o(yymJyJy.y;) 
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'^ui' 


+ 15(K,.Kim5KjK,r) + + 15(y,,Y^YjY* 

+ 9 o(K,.Kjmj;K„y;) + 9 o(K,,y,;y„m=y,;) + 9 o(K,.Kjny;'m^f,) + 8093 "Af; 
+ (y,.r,;) ( 49 ^^/! + sogiM, + 9 o(y,rT,.)) + i 5 (yJyiT;r,.) + 9 o(yJy. 

+ i5(Tjyjy;T,,) + 9o(riy„y,rT,,) + i5(rjyjy,,r,;) + 9o(Tjy„-y,,r;) 

- 49 f(r,,r,;) - 80 g|(r,.r,;)) 




ryn* rri 

uJ-1' 


UVo(2) 








— 545f2l^2p + 
64 


(B96) 

(B97) 


= + 32gti\Ms\^ - 2(2(f,t;)^(f;t,,)^ + 2(T;y*,)^ 

+ 2 (y:y,) (t,t;) + ^ ^ + 2 ^ ^ + 2 

) + 2(t;T*,)^ + 2(F;mf + (m^/Y;Y,, 

(B98) 


d-^t' 

J \ /I 


+ ) . (^2 




( 1 ) 


= —1 
25 


{2mgt\M:,\^ + 7gt\M^\^) 




aaS 


= 2mi„AY:Yt> 


= + 


ut' \ u 

256 


(B99) 

(BlOO) 


25 


gtim\^ + 32gll\Ms\^ 


2 

5 




mz^AYAYt, 


(30T„FJ(f*,T,i) +30F„Fj(F,,mjF;) +30 F.FJ(f*,^ 

+ 3oydTt + 30 F.FJ+ ^2^2,_ 15^2 

+ 15m2,, Jr(KFj) (fJF,,)^ + 20m^„ (f.F;)^ (f;F,,)^ + 10m? (f.F;)^ 

+ 5(f*,F;) (6(2m^^ + m?)F.y4 + 3(2T„Tt + 2F.m;Fj + + Y^Y^ml) 


+ 4m2,,^(F;F,,) ) + io(f„t,i)^(f;t*,) + io(f.f;)^(t:t,,) + io(t:f,,) ( 


TuYA 


+ 10 F:Ft 


T„T;) +5(Y:Yt, 


<y^ya 


+ 10(F;F,)^ (F„mjF;)^ + 10 (f;F,,)^ + 3{Y^YA)^[ml*Y:Y, 

+ io(KF;)^(y:mf F,,) + 5m|^-, ^(f;fJf;f*,) + 


y:yjy:y. 


^UV^(2) 
^ mf. 


{3Qgt\MA^ + 7gt\M^\- 


aa, 


( 2 ) 


= “l 
25 


(BlOl) 

(B102) 

(B103) 
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= (B104) 

- f .^IlMsr + (y.Y;) + 4m? (y,Y;) + 4(F,,mjF;) 

+ 2(F*,yJm2,,) + 4 (t,,T;) - 4^<7icri,i (B105) 

^3 ^ ^ ^ + iOglMs + SOglM,'^ - 90M, ) 

+ 5 (^Qgl ^ — 15(73M 3 + 8(7i ^ 2 M 3 + ?14i j j M| 

- 3 (240 + ml) ( (f,M;) f + 6gj [Y^m^Y;) - 90gl [Y^m^Y;) + 3gl (F,,Fjm|p) 

- A5gl [Y.Y^ml,,) + SOmj,^ {y,yJy,Y;) + 30m^„ {y,yJy,Y;) 

+ 30m| (^t'Y^YdY^) + 60m^^ + 30m| 

+ 3o(f*,T]T,F;) +3o(F,,TiT„F;) +3o(F,,mjF;F,F;) +3o(F,,mjFjF„F;) 

+ 30(F,,Fjm2F,F;) + 3o(F,,FjF,mjF;) + 15(F,M]F,Fjm2,,) 

+ 30(F*Mjm2F.F;) +30(F,MjF„mjF;) +15(F*,FjF„Fjm2p) +90(72 M*(f;T,,) 

+ 30(F;FdT;Tt,) +30(FjF„T;Ti,) +30(T]FrfF;Tt,) +30(TiF,F;Tt,) 

+ 3o(tJf;F,,T;) + 3o(tJf;F,,T;) + 6gl{T,,T:^ - 90gl{T,,T:^ - 80gta2,3 - 16 ( 73 ^ 2,11 

+ 8v^(7ia3,i + 90(F,,m5F;)Tr(F„Fj) + 45(Ft,Fjm2, )Tr(F,Fj) 

+ 90(Ti,T;)Tr(F„Fj) + (f,,T;) (24o(f;T,,) - Gg^M^ + 99glM2 + 90Tr(FjT„)) 

+ 90(Y;Te)Tr(T*Yl^) +6(Yt>Y;) {glm]j^ - IGglmjj^ + glmj - IGg^mj - 30^^|M2|^ 

+ 40(F*,m2F;)+30(F*,FX^-,) 

+ 40(t,,T;) + 30m^„Tr(F„Fj) + 15m?Tr(F„Fj) + 15Tt{t:T^) + 15Tr(mjFjF„) 

+ 15Tr(m2F„Fj))))) (B106) 


= 64(73^|M3p + (B107) 

+ 8gf\M,\^ - v^(7iai,i) (BIOS) 

/9i| = ^ (2(7^ (393(7 ?Mi + A9gl (2M1 + M3) ) Ml 

+ 5^2(73^ — IGg^M^ + 8g\{2M^ + Mi)) M 3 + 30(73(T2,3 + 3VT5(7i(T3,i + 6(73 (T 2 , 11 )) 

(B109) 
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(BllO) 


^mi~, . 

ut ,r 


= --(^(2glmH^ - SOglrriH^ + gimp - Ihglm}, + 4^i |Mip - + 80m^^ (YfYp^ 

+ 60m|(y,,y,;) + 4 o(y,m 2 y,r) + 4o(y,.y,;m|,-,) + 4o(r,,r,;) + 60m|,.Ti(y„y„t) 

+ i5m|TY(y„y„t) + 30Ti(r„Tj) + 30Ti(mJyiy„) + 30iv(mjy„y„t)) (y„-y,,)_ 

- 2jfAf;(y„T.,)_ + sogiM; (y.Tf), + 4o(y,.T;) (y„Tf)_ + 30Ti(r:y/) (y.Tc), 

- 2jfAf,(T;y,.)_ + 30glAk(T:Y,,'j^ + 4o(y,;T,,) (rjy.-), + 30Tr(y„tr„) (T;y,-)_ 

+ 2g;(r,:r,,)_ - 3092"(r,;r,,)_ + 4o(y,.i'r) (r„T,,)_ + 30Tr(i'„y„t) (r,:r,,)_ 

+ 9?(m?y„-y,-)_ - I5!,|(mry:y,.)_ + 2 o(y,,y;) (mry„-y,,)_ 

+ 15 iv(y.y„t) (mryjyg)^ + 2g't(Y;ml-Yg')^ - 3092 (y„-mJ-y,. 

+ 4o(y,.y;) (y,;mf y.), + 30Tr(y„yj) (y;mf y.), + lom 
+ lorni,, (y;y7y;y,.)_ + 5 m? (y;y/y;y)_ + io(y.;y/T,T,.)_ 

+ 20mi„ (y:yJy:y,,')_ + 5 m? (y;y.fy„-y,-)_ + io(y.;yjT.;r,.)_ 

+ io(y„TjT;y,,), + io(y;T,rr„*y,,)_ + io(r„-y/y;T,.)_ + io(r„-y„^y,:r,.)_ 

+ io(T„Tjyjy,,), + io(r:Tjy;y,,)_ + 5(m=*y,;y/yjy;.), + 5(m;-y„-yJ'y„-y,.)_ 

+ io(y„*mf y7y;y,,)_ + io(y„-mf yjy;;y.)_ + io(y;y7m5-y;y,,), + io(y„-y/y;m, 
+ io(y„*y„^mJ-y;y,,)_ + io(y„-y7y,;mf y.) j a‘'''/3® 


V* V 

''Hi \ 


= -Qgl\M2\^ + -i=gi(ypi - -^gl\Mi\^ - ^c/glMs] 


9 

^ 2| /I y |2 „2| /I y |2 




15^ 


(Bill) 

(B112) 


^ mi, 
Q 


— +gfi'lfi'2l^2p + 875f2l^2p + 325f25'3|M2p 


+ (^6(73 ^2Mi + M 3 j + 95^2 ^2Mi + M 2 j j + SGTgf^Mi j Mj" 

+ ^fi's (^13 ^ 105 f 3 M 3 + 3gl ^2M3 + M 2 j + gl j j "^3 + ~^9i92^i^2 

32 2 1 

+ l%2fi'3"^3^2 + 6fl'2‘^2,2 + “^fi'3<^2,3 + Y^fi'l<^2,ll + (B113) 


A 


= -6c/^|M2p - ^c/iCTyi - ^gllMil"^ - ^gllMsl'^ 


(B114) 



— +-^5'i5'2l^2p + 875f2l^2p + 325(2(73|M2p 

q' 0 

+ ^ 2 Mi + M 3 j + 9(72 (^2 Mi + M 2 j j + 867(7^Mi j M^ 

+ ~^ 9 i (^15 ^10(73M3 + 35'2 (^2M3 + M2 j j + (7i ^2M3 + ^1 j j ^3 + ~^ 9 i 92 ^ i ^2 


32 2 1 

+ 16(725'3^3^2 + 6(720'2,2 + -^5'3'^2,3 + Y^5'iC’'2,ii — 4 ^^^( 71(734 (B115) 

- 12^/i|Mi|2 + ^(Ti.i) (B116) 

^ ^5'i(l55fiO'2,ii + 4865(3|Mi|^ + 5A/l5cT3,ij (B117) 

/Si'], = -^^i(i2(7i|Mi|2 + (B118) 

/^S, = ^5'! (l%iO' 2 ,ii + 4865(3 |Mi|^ - 5 ^/ 150 - 3,1 j (B119) 


Vacuum expectation values 



= 3«.(y,.py,y;;) - |g;„. 


(B120) 

A‘'''/3l? = -/j(75g^ + 7gj)g. 


(B121) 

A/?« = -3g„(y,.y,;) 


(B122) 

(^5 (5 ( 32(72 + 9(72'Xi) + gl (oXi + 8) ) 

(nv,;) 

-45o((y;.y;))" 

+ 6( - 25(y,,y7riF;) + 2( - 75(y,,y;ty„y,;) 


) (B123) 



(B124) 
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